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INTRODUCTION 

The integrability of families of two degrees of freedom potentials with an invariant 
plane F = {x2 = 2/2 = 0} and normal variational equations (NVEs) of polynomial 
(degree two and odd degree) type along generic curves in T was studied jointly by 
the first and second author in [2 [5] . Related problems with this approach have 
been studied before by Morales and Simo (see [TT] ) and by Baider, Churchill and 
Rod (see [5]). The use of techniques of Differential Galois theory to determine 
the non-integrability of hamiltonian systems, appeared independently for first time 
in [71 [12] and [4], followed by |3], [1] and [11]. A common limitation presented 
in these works is that they only analyzed cases of fuchsian monodromy groups, 
avoiding cases of irregular singularities of linear differential equations. The case of 
the NVEs with irregular singularities can be approached from the Morales-Ramis 
framework ([Sl[n], see also [1]). 

In [2] , in the frame of Morales-Ramis theory, analyzing the NVEs through Ko- 
vacic algorithm, were obtained the following results: 

(1) Theorem. Galois group of the Hill-Schrodinger equation, 

with Pn{t) e C[t] non- constant polynomial, is a non-abelian connected group 
isomorphic either to SL{2,C) or to semidirect product of C* with C (also 
known as the Borel group). 
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(2) Non-integrability results. Let H be a hamiltonian system of two degrees 
of freedom given by H = T + V . If T is an invariant plane and the NVE 
around a generic integral curve in T is of Hill-Schrodinger type with degree 
2 polynomial coefficient, then the potential V is written such as follows: 

^ = T\ — V? + ^0 - ^1X2 - X2X1XI - X^x^xl + (3{xi,X2)xl, 

(A2 + 2X3X1)'' 

with A3 ^ 0; and thus they are not integrable by terms of rational functions. 
If the potential V is of the following form: 

V \o + Q{xi)xl + (3{xi,X2)xl, 

being Q{xi) a non-constant polynomial, then the it has NVE of Hill- 
Schrodinger type with non-constant polynomial coefficient along a generic 
integral curve in T — {x2 — y2 — 0}, and thus they are not integrable by 
terms of rational functions. 

Along this paper we follow [2] considering the NVEs of type Hill-Schrodingcr 
equation with cuartic polynomial potential, i.e cuartic polynomial normal varia- 
tional equations. We analyze deeply all the possibilities for V to be this kind of 
NVEs. 

1. MORALES-RAMIS THEORY 

Morales-Ramis theory relates the integrability of hamiltonian systems with the 
integrability of linear differential equations (see [8] [9] and see also ^). In such 
approach the linearization (variational equations) of hamiltonian systems along 
some known particular solution is studied. If the hamiltonian system is integrable, 
then we expect that the linearized equation has good properties in the sense of 
differential Galois theory (also known as Picard-Vessiot theory) . To be more precise, 
for integrable hamiltonian systems, the Galois group of the linearized equation 
must be virtually abelian. This gives us the best non-integrability criterion known 
so far for hamiltonian systems. This approach has been extended to higher order 
variational equations in |10] . 

1.1. Integrability of hamiltonian systems. A symplectic manifold (real or com- 
plex), M2n is a 271-dimcnsional manifold, provided with a non-degenerate closed 
2- form UJ2- This closed 2-form gives us a natural isomorphism between vector bun- 
dles, b: TM T*M. Given a function H on M, there is an unique vector field 
Xh such that, 

\>{Xh) = dH 

this is the hamiltonian vector field of H. Furthermore, it has a structure of Poisson 
algebra over the ring of differentiable functions of M2n by defining: 

{H,F} -.^XhF. 

We say that H and F are in involution if and only if {H, F} ~ 0. From our 
definition, it is obvious that is a first integral of Xh if and only if H and F are 
in involution. In particular H is always a first integral of Xh. Moreover, if H and 
F are in involution, then their flows commute. 
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The equations of the flow of Xh, hi a system of canonical coordinates, pi, . . . ,P' 
qi, . . . ,qn (that is, such that lu2 — J2^=i Pi ^ ^^^^ be written in the form 
fill f)M 
q = ^{={H,q}), P=-^h{H,p}), 

and they are known as Hamilton equations. 

Theorem 1.1 (Liouville- Arnold) . Let be a hamiltonian defined on a real sym- 
plectic manifold M2„. Assume that there are n functionally independent first inte- 
grals Fi, . . . , Fn in involution. Let Ma be a non-singular ( that is, dFi , . . . , dFn are 
independent over each point of Ma) level manifold, 

Ma = {p: Fi{p) = ai, . . . , F„(p) = a„}. 

(1) If Ma is compact and connected, then it is a torus Ma — K"/Z". 

(2) In a neighborhood of the torus Ma there are functions /i, .../„, 0i, 
such that 

n 

ti^2 = ^ dli A d4>i, 

i=l 

and {H, Ij} = for j — 1, . . . ,n. 

From now on, we will consider C^" as a complex symplectic manifold. Liouville- 
Arnold theorem gives us a notion of integrability for hamiltonian systems. A hamil- 
tonian H in C'^" is called integrable in the Liouville's sense if and only if there exists 
n independent first integrals of Xh in involution. We will say that H is integrable 
by rational functions if and only we can find a complete set of first integrals within 
the family of rational functions. 

1.2. Variational equations. We want to relate integrability of hamiltonian sys- 
tems with Picard-Vessiot theory. We deal with non-linear hamiltonian systems. 
But, given a hamiltonian H in C^" and T an integral curve of Xh, we can consider 
the first variational equation (VE), such as 

in which the linear equation is induced over the tangent bundle represents a 
vector field supported on F). 

Let F be parameterized by 7: 1 1— > {x(t), y{t)) in such way that 

dxi dH dyi dH 
dt dyi ' dt dxi 
Then the VE along F is the linear system, 

From the definition of Lie derivative, it follows that 

n rr 

dyi OXi 
is a solution of the VE. We can use a generalization of D'Alambert's method to 
reduce our VE (see [8l|9] and see also [6]), obtaining the so-called normal variational 
equation (the NVE). We can see that the NVE is a linear system of rank 2{n — 1). 
In the case of hamiltonian systems of 2-degrees of freedom, their NVE can be seen 
as second order linear homogeneous differential equation. 
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1.3. Non-integrability tools. Morales-Ramis theory is conformed by several re- 
sults relating the existence of first integrals of H with the Galois group of the 
variational equations (see for example [Sj, [S] and see also [B]). 

Most applications of Picard-Vessiot theory to integrability analysis are studied 
considering meromorphic functions, but for every equation considered throughout 
this paper, the point at infinity, of our particular solution, plays a transcendental 
role: the Galois group is mainly generated by Stokes phenomenon in the irregular 
singularity at infinity and by the exponential torus. So that we will only work with 
particular solutions in the context of meromorphic functions with certain properties 
of regularity near to the infinity point, that is, rational functions of the positions 
and momenta. Along this paper we will use the following result: 

Theorem 1.2 ([8 ). Let H he a hamiltonian in C^", and 7 a particular solution 
such that the NVE has irregular singularities at points of 7 at infinity. Then, if H 
is completely integrable by rational functions, then the identity component of Galois 
Group of the NVE is abelian. 

Remark 1.3. Here, the field of coefficients of the NVE is the field of meromorphic 
functions on 7. 

2. Method to determine families of hamiltonians with specific NVE 

This method was implemented in [2] as a generalization of the method shown in 
[llj . This section is devoted to this method. 

Let us consider a two degrees of freedom classical hamiltonian, 

V is the potential function^ and it is assumed to be analytical in some open subset 
of C^. The evolution of the system is determined by Hamilton equations: 

dV . dV 

0x1 0x2 

Let us assume that the plane F = {x2 = 0,2/2 — 0} is an invariant manifold of 
the hamiltonian. We keep in mind that the family of integral curves lying on F is 
parameterized by the energy h = H\r, but we do not need to use it explicitly. We 
are interested in studying the linear approximation of the system near F. Since F 
is an invariant manifold, we have 



-0, 

r 



dV 

dx2 

so that the NVE for a particular solution 

t ^ 7(t) = {xi{t),yi = xi{t),X2 = 0, y2 = 0), 

is written, 



Let us define. 



dxl 



{xi{t),0) 



8x2 
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and then we write the second order Taylor series in X2 for V ^ obtaining the fohowing 
expression for H 

(1) H=y^y^ + <j>{x^)-a{x^)'^+P{x^,X2)xl 

which is the general form of a classical analytic hamiltonian, with invariant plane 
r, providing that a Tayfor expansion of the potential around {x2 — 0} exist. The 
NVE associated to any integral curve lying on F is, 

(2) i^a{x,m- 

2.1. General Method. We are interested in computing hamiltonians of the fam- 
ily ([1]), such that its NVE ^ belongs to a specific family of Linear Differential 
Equations. Then we can apply our results about the integrability of this LDE, and 
Morales-Ramis theorem to obtain information about the non-integrability of such 
hamiltonians. 

From now on, we will write a{t) = a{xi{t)), for a generic curve 7 lying on F, 
parameterized by t. Then, the NVE is written 

(3) I = «we 

Problem. Consider a differential polynomial Q{ri, r],rj, . . .) e C[r], r), 77, . . be- 
ing 77 a differential indeterminate (Q is polynomial in rj and a finite number of the 
successive derivatives ofrj). Compute all hamiltonians in the family ^ verifying: 
for all any particular solution in T, the coefficient a{t) of the corresponding NVE 
is a differential zero of Q, in the sense that Q{a, d,a, . . .) — 0. 

We should notice that, for a generic integral curve 7(f) — {xi{t),yi = ii{t)) 
lying on F, ^ depends only of the values of functions a, and (f). It depends on 
a{xi), since a{t) — a{xi{t)). We observe that the curve j{t) is a solution of the 
restricted hamiltonian, 

(4) /i=| + </.(xi) 

whose associated hamiltonian vector field is, 

d dcp d 

(5) Xh=yi^ 1— 

0x1 axi oyi 

thus xi{t) is a solution of the differential equation, xi — —-j^, and then, the 
relation of xi{t) is given by cf). 

Since "f{t) is an integral curve of Xh, for any function defined in F we 

have 

l^*{f)^Y{X,f), 

where 7* denote the usual puUback of functions. Then, using a(t) = 7* (a), we have 
for each fc > 0, 

(6) ^=^*(^'")' 
so that, 

Q{a,a,a, . . .) = Q(7*(a), 7*(X/ja), 7*(X^a), . . .). 
There is an integral curve of the hamiltonian passing through each point of F, 
so that we have proven the following. 
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Proposition 2.1. Let H be a hamiltonian of the family ([T]), and Q(a, d,d, . . .) a 
differential polynomial with constants coefficients. Then, for each integral curve 
lying on T, the coefficient a{t) of the NVE ^ verifies Q{a, d,d, ...,) — 0, if and 
only if the function 

Q{xi,yi) = Q{a,Xha,Xla, . . .), 

vanishes on V . 

Remark 2.2. In fact, the NVE of a integral curve depends on the parameterization. 
Our criterion does not depend on any choice of parameterization of the integral 
curves. This is simple, the NVE corresponding to different parameterizations of 
the same integral curve are related by a translation of time t. We just observe that 
a polynomial Q{a, d,d, . . .) with constant coefficients is invariant of the group by 
translations of time. So that, if the coefficient a{t) of the NVE (0) for a certain 
parameterization of an integral curve ^(t) satisfied {Q — 0}, then it is also satisfied 
for any other right parameterization of the curve. 

Next, we will see that Q{xi,yi) is a polynomial in yi and its coefficients are 
differential polynomials in a, 0. If we write down the expressions for successive Lie 
derivatives of a, we obtain 

/r,N ,^ da 

(7) Xha = yi- — , 



(8) 



vl 



dxi 



d(j) da 
dxi dxi 



(9) 



,d^a 



dxj 



da 



dxi \dxi dxi 



d^a 



dxi dxf 



d^a 
dx\ 



vl 



d ( d ( d(f> da 
dx\ \dx\ \dx\ dx\ 



dx\ dx\ ) dxf dxi 



(10) 

In general form we have, 
(11) 





^ d4> da\ 


^ ^ d(i) d'^a\ 


i^dxi 


\dxi dxi J 


dxi dxi J 



dx- 



dxi dxi dyi 



it inductively follows that they all are polynomials in j/i, in which their coefficients 
are differential polynomials in a and cf). If we write it down explicitly. 



(12) 



n>k>0 



we can see that the coefScients En.k{ct, 0) G C 
recurrence law, 



satisfy the following 



(13) En+i.kia,4>) = 3 — £'„,fc_i(a,0) - (fc + l)£'„,fc+i(a, 0)-5 — 

dxi dxi 
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with initial conditions, 

(14) Ei,i{a, <!>) = —, ^i,fc(a,0) =0 Vfc^^ 1. 

axi 

Remark 2.3. The recurrence law ilS\) and the initial conditions Jj^p , determine 
the coefficients i?„^fc(a, 0). We can compute the value of some of them easily: 

• -B„,„(a,0) = ^ for alln>l. 

• -£'„.fc(a, (p) — Q if n — k is odd, or k < 0, or k > n. 

3. Main Result: Families of Hamiltonians with cuartic NVEs 

Theorem 3.1. Let H = T + V be a classical hamiltonian with invariant plane T 
such that the generic NVE along integral curves inV is a Hill- Schro dinger equation 
with cuartic polynomial coefficient. Then, the potential V is written in the form: 

V = Xo+ P4^{xi)xl + P{xi,X2)xl, 

where Aq G C, P4 is a cuartic polynomial and (3{xi,X2) is analytic around T. 

Proof. Following our general method, the family potentials of potentials satisfying 
the assumptions of the theorem are them given by solutions a{xi), (l){xi) of the 
system of differential equations 

£;5,5(a,0) =0, E5,3ia,(b)=0, ^5,1^ = 0. 

The first equation is just the following, 

d^a 

-^5,5 (a, 0) =0, 

and then we know that a is a cuartic polynomial in xi. 

a = a + bxi + cXi + dxi + exi- 

where 6, c, d, e are complex numbers and e does not vanish. Then we substitute a 
for his value in the equations 

£;5,3(«,0) =0, S5,l(«:0)=O, 

obtaining the following system of differential equations in 0(a;i): 

{Aexl + 3dxl + 2cxi + h)(j)''' + {mexl + SOdxi + 10c)<?i"'+ 

(L) + {24Qexi + 60d)(^" + 2406?^' = 0, 



{lM+72exi)((j)'f + {h + 2cxi+Mxl+Aexl){(i>"f+ 

(NL) + (14c + 'i2dxi + 84ea:^)0>" + (6 + 2cxi + Uxl + 4ex?)<?!)>"' = 0. 

Equations ([L| and (|NLp are ordinary differential equations in (j)' with some com- 
plex parameters. We can substitute a new unknown y for (j)' in order to reduce the 
order by one. Secondly, by a translation of xi by an scalar value, 

x — xi — 

we can assume that one of the coefficients of the polynomial ck(xi) vanish. From 
now on let us write, 

a{x) — a + bx + cx^ + ex^, 
and let us study the system of differential equations: 
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(L2) (4ex^ + 2cx + b)y"' + (GOex^ + 10c)y" + 240ea;?/' + 240e?/ = 0, 

(NL2) 72exy^ + {h + 2cx + Aex^){y'Y + (14c + Mex'^)yy' + (& + 2cx + Aex^)yy" = 0. 

The first equation (|L2[) is a linear equation in y. In this special case we will be 
able to completely solve the equation (jL2|) and then prove that solutions of (IL2|1 do 
not satisfy the non-linear equation ()NL2|) . Then, the only solution of the system is 
given by the function y = 0, that correspond to = Aq and then the potential 

V^^+^^xl+P{xi,X2)xl 

is of the form given in the statement of the theorem. 
Solution of the equation (|L2p 

This linear equation is solvable by elemental methods. Fortunately, its Galois 
group is trivial, and then we can look for a fundamental system of solutions that are 
rational functions over x, a, 5, c, e. The main problem, is that such a system of solu- 
tions do not always specializes to a particular system of solutions when determining 
the values of the parameters a, b, c, e. There are some values of the parameters that 
correspond to degeneration of the system of solutions. The equations of this locus 
of degeneration is given by the wronskian of the fundamental system. When the 
wronskian vanish, the fundamental system degenerates, and then there appear a 
different solution. We have to consider also these restricted problems independently. 



First, we find the general solution of l|L2|l . for generic values of the parameters, 
depending on arbitrary constants Ki, K2, K^: 

K,N, + K2N2 + KsN^ P 

(15) y = ^3 = 

where, 

D = Aex^ + 2cx + 6, 
Ni = xiiec^x^ - A2hecx^ - (6c^ -I- ^Seb'^)x^ + 9b^cx + 6b^) 
N2 = x{8ecx^ ~ 12bcex^ - (2466^ -I- 12c^)x^ - 12bc^x^ + 36^) 
N3 = Sc^e'^x^ - Mbce^x^ - {I2c^e + im'^e^)x^ + 21b^ex - 3b^c^ 

Let us study for which parameters the above expression is not the general solution 
of (|L2[) . It happens if and only if the wronskian of the fundamental system of solu- 
tions vanish. We know that the wronskian of the fundamental solutions Ni / D van- 
ish if and only if the wronskian of the numerators Ni vanish. We compute it, obtain- 
ing l^(iVi, TVa, iVs) = 162c367 + i2966Va;-H 38886^0^x2 -|-(259266ec3 -1-51846^)2:3-1- 
(2592c^63 -t- 1555265c4e)a;'' + ■^llQAb^e'ex^ + {1^^^2b^c'e^ + 2Q7mb^c^e)x^ + 
62208b^c'^e^x'' + 62208b^c'e^x^ + 414725V e^a;^ + 829 Ub^c'^e^x^° + AU72b^c^e*x^^ . 
The equations of the vanishing of the wronskian are given by the coefficients of this 
polynomial. This system is simple to solve and it has two independent solutions, 
that we will consider independently. 

(a) {b = 0} 

(b) {c = 0} 
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Case A, 6 = 0. 

If b vanish then the considered system of equations is: 
(L3) (4ea;3 + 2cx)y"' + {60ex^ + 10c)y" + 2A0exy' + 2AQey = 0, 

(NL3) 72exy^ + (2cx + Aex'^){y'f + (14c + Mex^)yy' + {2cx + Aex^)yy" = 0. 

we obtain a new general solution for this restricted case by direct integration. The 
general solution is: 

(16) y^ 

where, 



D3 = 2ex^ + c, 
^3,1 — 6ex^ — c, 
iVg 2 = x{-3c + 2ex'^) 
3 ^ (c3 + 6ec^x^ + 16e^x^)x-^ 
In this case, the wronskian of the numerators is, 

, 96ec(6ecV + IGe^x^ + 5c3) 

W 7V3,i,iV3,2,iV3,3 = ^ 4 ^, 

x^ 

so that this system of solutions only degenerates when c = which is considered as 
a particular case of the following. 

Case B, c — 0. 

Let us finally consider the last case c — 0. The system of equations is now as 
follows. 

(L4) (4ea;3 + b)y"' + GOex^y" + 2A0exy' + 2A0ey = 0, 

(NL4) 72exy^ + (6 + Aex^){y'f + SAex^yy' + {b + 'iex^)yy" = 0. 

The general solution of (|L4p is given by, 

^''^ ' = 151 " = ^ 

where 

D4 = iex^ + b, 
Ni^i = x{b - 8ex^), 
Ni^2 = x^{b~2ex^) 
Ni,3 = [h^ - 2Sebx^ + lee^x^). 
We analyze the vanishing of the wronskian, Vl^(iV4,i, iV4,2, -^4,3) = 2b'^ + 32eb^x^ + 
192e'^b^x'^ + 512be^x^ + bUe^x^'^. Because of the coefladent in the 12th power of X, 
this wronskian does not vanish for any of the considered values of the parameters, 
and then this general solution does not degenerates. 

Common solutions with the non-linear equation 

Here we look for solutions of the linear equation that also satisfy the considered 
non-linear equation. We directly substitute the general solution (jisp of (|L2[) into 
(|NL2[) . Then we obtain an rational expression: 

Q{x;b, c,e,Ki,K2,K3) _ 
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being Q a polynomial in x depending on the parameter b, c, e, Ki, K2, K3. Thus, 
we look for the values of the parameters that force Q to vanish. If we develop Q 
as differential polynomial in P and D we obtain the following expression: Q = 
{72exD + (-6 + 42c + 262ex^)DD' + 21D'^ + ~ 3DD")P'^ + (14c + 8Aex)D^ - 
6DD')PP' + D^P'^. Note that the polynomials D,D',D" does not depend on 
the parameters Ki, and P, P' are linear in such parameters. It follows that Q is 
polynomial in x of degree 16 whose coefficients are homogeneous polynomials of 
degree two in the parameters Ki, K2, K3. 

16 

Q = ^ C,(ifi, X2, if3, 6, c, = 0, 

1=0 

where: 

/All X\,\ (kA 

a = (ifl,if2,i^3) A^21 AI2 AI3 X2 . 

V\i AI2 \h) {KbJ 
Finally, the coefficients A^j. are polynomials in the parameters b, c, e. The common 
solution of (|L2p and (jNL2[) correspond to values of the parameters 6, c, e, Ki, K2, K3 
that are solutions of the system of 17 algebraic equations: 

aiKi,K2,K3,b,c,e)^0, z = 0,...,16. 

Each equation = is the equation of a cone 3-dimensional the affine space 
over the field C(6, c, e); thus the equation of a conic curve in the projective plane 
P(C(5, c, e); of homogeneous coordinates Ki, K2, K3. Two conic curves intersect in 
4 points. This simplifies the computations, that have to be carried symbolically, 
and lead to the incompatibility of the system, provided that b and c are different 
from zero. 

Following the same schema, we analyze the exceptional case 6 = 0. We put the 
solution of (|L3|) into the equation (jNL3p . We obtain an expression, 

x^Dl 

where, Q3 = {72exD3 + (-6 + 42c + 2Q2ex^)D3D'^ + 21D^2 + 1)2 _ SD3D'^)Pi + 
(14c + 84ex)Dl - 6D3D'3)P3P^ + D^P^. Here, we find that Q3 is a polynomial in 
X of degree 18, 

18 

03 = 5I^*(c,e,i^i,i^2,if3K, 

i=Q 

and again the system of algebraic equations {Ei = 0}, is qualitatively similar to the 
above system {Ci = 0}. The same analysis is carried out in the other exceptional 
case c = 0. In this last case we obtain that Q4 = {J2exDi + (—6 + 2(!)2ex^)D4^D'^ + 
21D2 + Dl- iDiD'DPl + 84ea;D| - &DiD'^)PiP!^ + n|P| is of degree 18 in x, and 
then we have a system {Fi = 0} 19 algebraic equations that form an incompatible 
system. 

□ 

Final comments and open questions 

One open problem presented in [2 is the problem of determining families of clas- 
sical hamiltonians with an invariant plane and NVE of Hill-Schrodinger type with 
polynomial coefficient of even degree greater than two, now the case of even degree 
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greater than four is still open. 

The problem of analyzing the monodromy of the NVE of integral curves of a 
two degrees of freedom hamiltonian (both, classical and general) has been studied 
by Baider, Churchill and Rod at the beginning of the 90's (see [3]). Their method 
is quite different, they have imposed the monodromy group to verify some special 
properties that were translated as algebraic conditions in the hamiltonian functions. 
Their theory were restricted to the case of fuchsian groups, which in terms of 
Galois theory means regular singularities, while we work in the general case. The 
comparison of both methods should be done. 

Another problem is to apply these method to higher variational equations. 
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